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摘 要 : 本 文 研究 了 一 类 非 线性 时 滞 系 统 的 Hoo 故障 检测 观测 器 设计 问题 。 以 状态 观测 器 作为 残 差 产生 
器 ， 建 立 观测 误差 动态 系统 。 在 系统 的 初始 状态 为 不 确定 但 有 界 的 情形 下 ， 对 这 类 观测 器 的 设计 
提出 两 类 改进 的 Hoo 性 能 指标 ， 并 以 线性 矩阵 不 等 式 (LMI) 形式 给 出 了 故障 检测 观测 器 的 设计 


方法 。 
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1 引言 


随 着 生产 过 程 的 日 益 大 型 化 和 复杂 化 ， 提 高 控制 系统 的 可 靠 性 变 得 尤为 重要 。 近 年 来 ， 
动态 系统 的 故障 检测 研究 引起 了 相当 多 的 关注 ， 一 些 不 同 的 故障 检测 方法 相继 被 提出 ， 如 神 
经 网 络 方法 申 ， 奇 偶 方程 方法 由， 基于 观测 器 的 方法 B"9 等。 事实 上 ， 大 多 数 实 际 控制 系统 
都 存在 非 线 性 性 ， 同 时 被 控 对 象 可 能 会 出 现时 滞 、 不 确定 性 或 外 界 干扰 。Takagi 和 Sugeno 提 
出 了 一 类 由 IF-THEN 规则 描述 的 模糊 动态 模型 ， 它 能 够 逼近 一 大 类 非 线性 系统 站。 近 十 年 
K., T-S 模糊 系 统 已 成 为 非 线性 控制 领域 中 的 一 个 重要 研究 方向 &9]。 最 近 ， 文 献 [6] 研究 了 不 
确定 T-S 模糊 系统 的 故障 检测 问题 ， 其 中 假设 对 象 初始 状态 为 零 ， 没 有 考虑 初始 状态 对 故障 检 
测 的 影响 。 一 般 地 ， 实 际 控制 系统 的 初始 状态 是 未 知 的 。 在 这 种 情形 下 ， 本 文 针对 由 工 -$ 模糊 
模型 表示 的 非 线性 时 濡 系统， 提出 由 故障 到 残 差 和 由 干扰 到 残 差 的 两 个 改进 的 五。。 性 能 指标 ， 
并 以 LMI 形 式 给 出 了 故障 检测 观测 器 的 设计 方法 。 


2 系统 描述 
考虑 如 下 一 类 由 T-S 模糊 模型 描述 的 带 有 于 扰 和 故障 的 非 线 性 时 滞 系 统 
z(t) = È hi(0(t)) [Aizx(t) + Aniz(t — o) + Biw(t) + Gif(t)], 
y(t) = S hi(O(t)) [Cia (t) + Cuiz(t — o) + Diw(t) + If), (1) 
a(t) = ol), t € [-o0,0], 
其 中 
hi(O(t)) 20, ¢=1,2,---,k, SoA, (O(t)) =1. 
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o > 0 表示 不 确定 的 定常 时 洁 ， 满 足 0 < o < co， 其 中 oo 为 已 知 常数 。z(t) € R” 和 yb e 
R! 分 别 是 状态 向 量 和 测量 输出 向 量 。w(t) e R? 和 f(t) e R 分 别 是 系统 的 外 界 干扰 向 量 和 故 
E, H w(t), f(t) € L2[0, 00)。 Ai, Au, Bi, Ci, Cii, Di, G 和 Ji (i = 1,2,..- k) 是 具有 适 
当 维 数 的 已 知 常 矩 阵 。 

实际 上 ， 系 统 (1) 的 初始 状态 一 般 为 不 确定 的 ， 但 属于 某 个 有 界 范围 。 因 此 ， 本 文 假设 : FF 
在 已 知 定常 矩阵 Mı, Mo E R”*?”, AM > 0, Mə >0, 对 任意 的 t € [~o0, 0], 有 


y(t) E€ ù = {2:27Myz<1,zER"}, plt)eN2={z2:z M2z<1, zER”}. (2) 
针对 系统 (1)， 本 文采 用 如 下 的 % 阶 模糊 观测 器 作为 残 差 产 生 器 
é(t) = 3 COES E A OE O 


k 
H(t) = D milO(e) [C2(t) + Cult- 0)], (3) 


a(t)=0, te [-o0,0], 


其 中 2&(t) EMMBRASAE, g(t) Eye 的 估计 ，r(t) = y(t) -9(t) 是 残 差 信号 。 观 测 器 增 
益 Lj (i =1,2,--- ,kk) 是 待 设计 的 。 简 记 


z(t— 0o) = zo(t), L(t—0o)=$(t), hi(0(t)) = hi. 


观测 误差 为 
e(t)=2(t)—#(), €(t) = zo (t) — ĉo (t), 


由 式 (1), (3) 得 观测 误差 动态 系统 为 


é(t) = = 5 hih; {(Ai 一 L;C;)e(t) + (Aii = LiC1;)eo(t) 
+ (Bi — LiD w(t) + (Gi — Lids) f(0}, a 


elt) = y(t), te [00,0]. 
AMR RZ S r(t) 的 Ls 范 数 作为 残 差 估 计 函 数 ， 即 
Ta 
f rT(t)r(t)dt, 
0 


EP Ty —TARIER. MET RAT RRM, MAO in. RE MARAE 
障 (f(t) = 0) 运行 状态 下 ， 残 差 估计 函数 的 最 大 值 。 常 用 的 故障 检测 的 方案 是 


Ta 
I r7 (t)r(t)dt > Jin > 有 故障 > 产生 警报 ， 
0 


Ta 
| rT (t)r(t)dt < Jon > 无 故障 . 
0 
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3 Hy 故障 检测 观测 器 的 设计 


引 理 10ol 对 于 任意 定常 矩阵 M > 0， 标 量 > 0， 和 向 量 函 数 z(t) : [0,0] > R” 在 下 述 
积分 有 意义 的 情况 下 ， 有 下 列 不 等 式 成 立 


Lf. z(s)as] "M | 2(s)ds| <o is zT (s)Mz(s)ds. (5) 
3.1 无 故障 和 有 干扰 情形 时 的 五 观测 器 


由 于 系统 的 不 确定 初始 状态 对 残 差 有 一 定 影响 ， 因 此 ， 本 文 将 设计 一 个 形 如 式 (3) 的 观测 
器 ， 满 足 性 质 














Ta Ta 
1 rT(t)r(t)dt < 7? J wT (t)w(t)dt + a1, (6) 
0 0 
其 中 和 ol 是 给 定 的 正常 数 。 考 


Ta 
J wT (t)w(t)dt < N, 
0 


NN 为 一 个 常数 ， 则 故障 检测 的 阔 值 J 选 为 Jein = YN +016 

注 1 与 文献 [0 中 设计 故障 检测 系统 的 瑟 w 指标 不 同 ， 式 (6) 为 一 个 改进 的 及 指标 ， 其 
中 ?2 表示 干扰 对 残 差 的 影响 指标 ，aa 表示 系统 初始 状态 对 残 差 的 影响 指标 。 

定理 1 考虑 无 故障 和 有 干扰 的 系统 (1)， 给 定 标量 由 ait, 7 > 0， 如 果 存 在 和 矩阵 已 > 
0, Qi (i =1,2,--- k), T, >0, To >0, 标量 041 >0, ai2 > 0 和 ala3 > 0， 使 得 LMIs 


Aij < 0， ij = 1,2,... ,k, (7) 
P-—ayM, < 0 T- ayoM, <0, To — @iz Mə <0, (8) 
Q11 十 aoaaz + of013 < a1, (9) 


Pee a FLO a 23 38 ah EE AL, = P71Qi(i=1,2,… ,k)， 其 
Aij 2 


Si — Vij +Ti + Qeiey PAn — QiCig+Qecierg dATP —dCP QT 0 PB; — QiD; + Noi,dy 
* Ne1i,e1g — Ti dAT, P — acT QF 0 Qe1i,dj 
* * ooT2 — 2dP 0 dP B; 一 dQiD; 
1 
* * * 一 So 12 0 
$ * * * Qai,d; — VT 


且 
Zi = PA + ATP, Vij =QiCj+CT QT, Reij = (CPC) +C7TCi)/ 2， 
Neier; = (CF C1; + CF Cu)/ 2, Qeri,erg = (CC + CYC1i)/ 2, 
Neij = (CTD; +CT D:)/2 Nosg = (CHD; + CF D:)/2, Naia = (DTD; + DT D,)/2. 


证 明 选择 Lyapunoyv 函数 


V(e(t)) = e? (t)Pe(t) + J eT (s)Tie(s)ds + T (s — t+ oo)ėT(s)Toė(s)ds, (10) 


t~o 
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其 路 和 TD 是 正定 矩阵 。V(e(t)) 沿 误差 系统 (4) 的 轨迹 对 时 间 求 导 得 
V (elt)) = 2eT(t) Pe(t) + eT (Telt) — eT(t)Ty eolt) 


+oé? (t)T2é(t) 一 a é7 (s)Toé(s)ds, 


进一步 ， 有 
V (e(t)) +r (t)r(t) — r7 (t)r(t) 


k k 
= P P hih fe" (0 [PAs — L0) + (AP - OP LPP +N pee Gs 


i=1 j=1 


CC i CC 


CF Cu + CP Cu 
En e(t) + eT) (AR — OL )P + “Je(t) 


+e? (t)| P(A1i — LiCij) + 


DT Cy 十 DY | 
2 


pe Et De) + wr [BF — DIL)P+ e(t) 


D7 D; + Dj Di 
Aio) 


CED; TELD: ‘Jul OW + ura [ZEP Oer] 


[ 

+e" (t) [PUB -L = 

tes | 
eae he 


十 ez (t) CHiCy > CHORE T|eo(t) + aoé? (t)Tré(t) — f i è (s)Toè(s)ds} —rT(t)r(t). (11) 


对 任意 正定 矩阵 已 > 0， 有 下 列 等 式 成 立 
k k 
0 = 27P, 》 > hihj{(4i — LiC;)e(t) 


j>] j=1 
+(A1; — LiCi;)ec (t)+(B; — LiD;)w(t)} — 2é7 Pre. (12) 
将 (12) 加 到 (11) 上 ， 应 用 引 理 1， 令 已 =dP, PLi = Qi 和 


a(t) = [er paoro fe | 


V (e(t)) < 人 5 hghjé? (t)Aijé(t) — rT (t)r(t) + yw" (t)w(t). 


i=l 7 一 1 


由 条 件 (7)， 有 
V (e(t)) < rT (t)r(t) + wT (t)w(t), 


对 该 式 两 边 积 分 ， 有 
V(e(Ta)) - V(e(0)) < 矿 {=r (t)r(t) + uT (u(t) }dt. 
由 于 对 任意 的 te [-00,0], He(t) = z(t), elt) € Q, élt) € Q2， 则 由 (8), (9)， 得 
V(e(0)) = e7(0)Pe(0) + A : eT (t)Tye(t)dt + 人 (t + o0)é7 (t)T2é(t)dt 


2 
< a11 + 09012 + 99013 < Q1. 
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由 于 V(e(T4)) > 0， 得 
Ta Ta 
J r? (t)r(t)dt < vf wT (t)w(t)dt + a1, 
0 0 


即 所 设计 的 观测 器 满足 (6)。 
3.2 ”有 故障 和 无 干扰 情形 时 的 1 观测 器 
类 似 3.1 节 ， 针 对 对 象 的 不 确定 初始 状态 ， 设 计 一 个 形 如 (3) 的 故障 检测 观测 器 ， 满 足 


Ta Ta 
i rT(t)r(t)dt > B? f FOr — a. (13) 
0 0 


同样 地 ，(13) 也 是 一 个 改进 的 五, 指标， 其 中 8 和 as 是 给 定 的 常数 。B? 表示 故障 对 残 差 的 
影响 指标 ，as 表示 对 象 初 始 状态 对 残 差 的 影响 指标 。 

定理 2 考虑 有 故障 和 无 干扰 的 系统 (1)， 给 定 标 量 d, az, 6 > 0， 如 果 存 在 矩阵 已 > 
0，Qi (i = 二 1,2,… ,kk), Tr > 0, Te > 0， 标量 al > 0,Q22 > 0M az > 0， 使 得 LMIs 


Ai < 0, i j =1,2,--- ,k, (14) 
P—aqM, <0, Ti —a92M, <0, Dh- &œz23Mə < 0, (15) 
Q21 + C0Q22 + ohang < QQ (16) 


成 立 ， 则 观测 器 (3) 满足 (13)， 且 观测 器 增益 矩阵 为 L; = P-1Q; (i =1,2,---,k)) HAAG A 


Ei — Vij + Tr — doc PAn — QiC1j — Qeiergy dATP — ACP QT 0 PG; — Qid; — Qei, 35 
* 一 Qu ci -Ti dAT,P — dC; QT 0 15,55 
* * ooT2 — 2dP 0 dPG; — dQiJd; 
* * * mee 0 
* * * * 74,55 + BPI 


A 
Noa; = (CPI; + CTH) /2, Reniss = (CRIs + CU) /2, Quisi = (JEJ; + JF Ji) /2, 


Bi, Wiz, Duo Reij 和 foi,c1; 同 定理 1。 

证 明 “类似 于 定理 1 的 证 明 ， 证 略 。 

3.3” 既 有 干扰 又 有 故障 情形 时 的 Hoo 观测 器 

定理 3 ”考虑 既 有 干扰 又 有 故障 的 系统 (1), 给 定 标 量 d, a, co, 如 果 存 在 矩阵 P > 0, Qi (i = 
1,2,---,k), T > 0, T2>0, 标量 Qs1 > 0, age > 0l azz > 0， 使 得 以 下 的 优化 问题 


Min : J(7,8) = Y-B, (17) 
(7), (14), 

s.t. LMIs : P+ &zı Mı < 0, Ty = 32M, < 0, 了 2 = &33 Mə < 0, 
a31 + 79032 + 06033 < a, 


有 解 ， 则 观测 器 增益 矩阵 为 L; = P-1Q; (i = 1,2,… ,k)， 且 所 设计 的 观测 器 满足 (6) 和 (13), 
即 观测 器 对 外 部 干扰 的 抑制 度 为 72， 同 时 对 故障 的 灵敏 度 为 82。 
证 明 令 ai = as =aw， 由 定理 3.1 和 定理 3.2， 以 及 优化 问题 的 目标 函数 ， 易 证 。 
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Fault Detection for a Class of Nonlinear Delay Systems 
—An LMI Approach 


ZHANG Yu, JIA Xin-chun, ZHENG Li-hong, ZHANG Da-wei 


(School of Mathematical Science, Shanxi University, Taiyuan 030006) 


Abstract: The design problem of the Hə fault detection observer for a class of nonlinear delay 
systems is studied. A state observer is considered as the residual generator and the dynamic system of 
the observation error is constructed. Two improved Hæ performance indexes for such observer design 
are proposed under the assumption that the initial states of the system are uncertain but bounded. 
Design methods for fault detection observer are given in terms of linear matrix inequalities. 
Keywords: T-S fuzzy system; delay; fault detection; fuzzy observer; linear matrix inequality 
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